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In this article, we investigate the structures of the pseudoscalar mesons (it, K, D, D s , B 
and B s ) in the framework of the coupled rainbow Schwinger-Dyson equation and ladder Bethe- 
Salpeter equation with the confining effective potential (infrared modified flat bottom potential). 
The Schwinger-Dyson functions for the u, d and s quarks are greatly renormalized at small momen- 
tum region and the curves are steep at about q 2 = IGeV 2 which indicates an explicitly dynamical 
symmetry breaking. The Euclidean time fourier transformed quark propagators have no mass poles 
in the time-like region which naturally implements confinement. As for the c and b quarks, the 

■ current masses are very large, the renormalization are more tender, however, mass poles in the 
time-like region are also absent. The Bethe-Salpeter wavefunctions for those mesons have the same 
type (Gaussian type) momentum dependence and center around small momentum which indicate 

£N| ' that the bound states exist in the infrared region. The decay constants for those pseudoscalar 

^ \ mesons are compatible with the values of experimental extractions and theoretical calculations, 

ON ■ such as lattice simulations and QCD sum rules. 
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I. INTRODUCTION 

The low energy nonperturbative properties of quantum chromodynamics (QCD) put forward a great challenge 
for physicists as the strong SU(3) gauge coupling at small momentum region destroys the perturbative expansion 
approach, which is popular in theoretical physics. The physicists propose many original approaches to deal with 
the long distance properties of QCD, such as Chiral perturbation theory [1], heavy quark effective theory [2], QCD 
sum rules [3], lattice QCD [4], perturbative QCD [5], coupled Schwinger-Dyson equation (SDE) and Bethe-Salpeter 
equation (BSE) method [6], etc. All of those approaches have both outstanding advantages and obvious shortcomings 
in one or other ways. The coupled rainbow SDE and ladder BSE have given a lot of successful descriptions of the 
long distance properties of the low energy QCD and the QCD vacuum (for example, Refs. [7-10] , for recent reviews 
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one can see Refs. [11,12]). The SDE can naturally embody the dynamical symmetry breaking and confinement which 
are two crucial features of QCD, although they correspond to two very different energy scales [13,14]. On the other 
hand, the BSE is a conventional approach in dealing with the two body relativistic bound state problems [15]. From 
the solutions of the BSE, we can obtain useful information about the under-structure of the mesons and obtain 
powerful tests for the quark theory. However, the obviously drawback may be the model dependent kernels for the 
gluon two point Green's function and the truncations for the coupled divergent SDE and BSE series in one or the 
other ways [21]. Many analytical and numerical calculations indicate that the coupled rainbow SDE and ladder 
BSE with phenomenological potential models can give model independent results and satisfactory values [6-12]. The 
usually used effective potential models are confining Dirac S function potential, Gaussian distribution potential and 
flat bottom potential (FBP) [11,12,16-18]. The FBP is a sum of Yukawa potentials, which not only satisfies gauge 
invariance, chiral invariance and fully relativistic covariance, but also suppresses the singular point which the Yukawa 
potential has. It works well in understanding the dynamical chiral symmetry braking, confinement and the QCD 
vacuum as well as the meson structures, such as electromagnetic form factors, radius, decay constants [19-21]. 

The decay constants of the pseudoscalar mesons play an important role in modern physics with the assumption 
of current-meson duality. The precise knowledge of the those values (especially the Jd , Jd s , fs and /b s ) will provide 
great improvements in our understanding of various processes convolving the D, D s , B and B s mesons decays. At 
present, it is a great challenge to extract the values of the B and B s mesons decay constants Jb and fs s from the 
experimental data while the experimental values for the D and D s mesons decay constants fo and Jd 3 are plagued 
with large uncertainties. So it is interesting to combine those successful potentials within the framework of coupled 
SDE and BSE to calculate the decay constants of the pseudoscalar mesons such as 7r, K, D, D s , B, and B s while we 
prefer the detailed studies of the quarkonium to another article. For previous studies about the electroweak decays 
of the pseudoscalar mesons with the SDE and BSE, one can consult Refs. [6-12]. In this article, we use an infrared 
modified flat-bottom potential (IMFBP) which takes the advantages of both the Gaussian distribution potential and 
the FBP to calculate the decay constants of those pseudoscalar mesons. The present article is a continuation of our 
previous work [21]. 

The article is arranged as follows: we introduce the infrared modified flat bottom potential in section II; in 
section III, IV and V, we solve the rainbow Schwinger-Dyson equation and ladder Bethe-Salpeter equation, explore 
the analyticity of the quark propagators, investigate the dynamical symmetry breaking and confinement, finally obtain 
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the decay constants for those pseudoscalar mesons; section VI is reserved for conclusion. 

II. INFRARED MODIFIED FLAT BOTTOM POTENTIAL 



The present techniques in QCD calculation can not give satisfactory large r behavior for the gluon two point 
Green's function to implement the linear potential confinement mechanism in practical manipulation, the phe- 
nomenological effective potential models always do the work. As in our previous work [21], we use a gaussian 
distribution function to represent the infrared behavior of the gluon two point Green's function, 



which determines the quark-antiquark interaction through a strength parameter w and a ranger parameter A. This 
form is inspired by the 5 function potential (in other words the infrared dominated potential) used in Refs. [16,17], 
which it approaches in the limit A — > 0. For the intermediate momentum, we take the FBP as the best approximation 
and neglect the large momentum contributions from the perturbative QCD calculations as the coupling constant at 
high energy is very small. The FBP is a sum of Yukawa potentials which is an analogy to the exchange of a series of 
particles and ghosts with different masses (Euclidean Form), 



where N stands for the minimum value of the masses, p is their mass difference, and aj is their relative coupling 
constant. Due to the particular condition wc take for the FBP, there is no divergence in solving the SDE. In its three 
dimensional form, the FBP takes the following form: 



4 7 rG 1 (fc 2 ) = 37r 2 ^e-^ 



(1) 



n 




(2) 




(3) 



In order to suppress the singular point at r — 0, we take the following conditions: 



V(0) = constant, 



dV(0) _ f_ 
dr i 



2 V{0) 
dr 2 



d n V(0) 



dr 



= 0. 



(4) 



So wc can determine a,j by solve the following equations, inferred from the flat bottom condition Eq.(4), 



Here we correct a writing error in our article [21] i.e. change "small r" to "large r 
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n 

J2a j (N + jp) = V(0), 

3=0 

n 

5>(iV + . 7 p) 2 = 0, 
1=0 

n 

Y / a j (N + j P r = 0. (5) 

3=0 

As in previous literature [18-21], n is set to be 9. The phenomenological effective potential (infrared modified flat 
bottom potential) can be taken as G(k 2 ) — G\(k 2 ) + G2{k 2 ). 

III. SCHWINGER-DYSON EQUATION 

The Schwinger-Dyson equation can provide a natural framework for investigating the nonperturbative properties 
of the quark and gluon Green's functions. By studying the evolution behavior and analytic structure of the dressed 
quark propagators, we can obtain valuable information about the dynamical symmetry breaking phenomenon and 
confinement. In the following, we write down the SDE for the quark propagator, 

1 (\tt f d 4 k 

S-\p) = i 1 .p + m+ — J —T^SikhvG^k-p), (6) 

where 

S-'ip) = zA(p 2 ) 7 ■ p + B{p 2 ) ee A(p 2 )[i 7 ■ p + m(p 2 )], (7) 
G„„« - (<W - ^)G(fc 2 ), (8) 

and m stands for an explicit quark mass-breaking term. In this article, we take the rainbow approximation =7^. 
With the explicit small mass term for the u and d quarks (comparing with the u and d quarks, the mass term for 
the s, c and b quarks is large), we can preclude the zero solution for the B(p) and in fact there indeed exists a small 
bare current quark mass. In this article, we take Landau gauge. This dressing comprises the notation of constituent 
quark by providing a mass m(p 2 ) = B(p 2 )/A(p 2 ), for the u, d and s quarks, which is corresponding to the dynamical 
symmetry breaking. Because the form of the gluon propagator G(p) in the infrared region can not be exactly inferred 
from the SU(3) color gauge theory, one often uses model dependent forms as input parameters in the previous studies 
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of the rainbow SDE [6,11,12,18-22], in this article we use the infrared modified FBP to substitute for the gluon 
propagator. 

In this article, we assume that a Wick rotation to Euclidean variables is allowed, and perform a rotation analyti- 
cally continuing p and k into the Euclidean region where they can be denoted by p and k, respectively. Alternatively, 
one can derive the SDE from the Euclidean path-integral formulation of the theory, thus avoiding possible difficulties 
in performing the Wick rotation [23] . As far as only numerical results are concerned, the two procedures are equal. 
In fact, the analytical structures of quark propagators have interesting information about confinement, we will make 
detailed discussion about the u, d, s, c and b quarks propagators respectively in section V. 

The Euclidean rainbow SDE can be projected into two coupled integral equations for A(p 2 ) and B(p 2 ), the 
explicit expressions for those equations can be found in Ref. [19,20]. For simplicity, we ignore the bar on p and k in 
the following notations. 

IV. BETHE-SALPETER EQUATION 

The BSE is a conventional approach in dealing with the two body relativistic bound state problems [15]. The 
quark theory of the mesons suppose that the mesons are quark and antiquark bound states. The precise knowledge 
about the quark structures of the mesons will result in better understanding of their properties. In the following, we 
write down the ladder BSE for the pseudoscalar mesons, 

1 PlTT f H 4 k 

S^( q + CP)x(q,P)SZ 1 (q-(l-m = —J j-y A l» X {k,P)l v G»v{q-k), (9) 

where S(q) is the quark propagator, G^ u (k) is the gluon propagator, is the four- momentum of the center of mass 
of the pseudoscalar mesons, is the relative four-momentum between the quark and antiquark in the pseudoscalar 
mesons, 7 M is the bare vertex of quark-gluon, £ is the center of mass parameter which can be chosen to between and 
1, and \(q, P) is the Bethe-Salpeter wavefunction (BSW) of the bound state. 

We can perform the Wick rotation analytically and continue q and k into the Euclidean region § , the Euclidean 

pseudoscalar BSW x(Q> P) can be expanded in lorentz-invariant functions and 5*0(4) eigenfunctions, Tchcbychev 

i 

polynomials T„ 2 (cos 6). In the lowest order approximation, the BSW x(g, P) takes the following form, 



s To avoid possible difficulties in performing the Wick rotation, one can derive the BSE from the Euclidean path-integral 
formulation of the theory. As far as only numerical results are concerned, the two procedures are equal. 
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X(«, P) = 75 [»*?(«, P) + 7 • ^ 2 °fe n + 7 • OT • PF£(q, P)+i[ 1 -q, 1 - P]P 4 °(g, P)] . (10) 

It is important to translate the wavefunctions P™ into the same mass dimension to facilitate the calculations in solving 
the BSE, 

F° - F?, P 2 ° - A 1 ^ , F3 1 - A 3 ^ 1 , P 4 ° - A 2 P 4 °, g - g/A, P ^ P/A , (11) 

where A is some quantity of the dimension of mass. Then the ladder BSE can be projected into the following four 
coupled integral equations, 

POO P7T 

H(l,l)F?(q,P) + H(l,2)F°(q,P) + H(l,3)Fi(q,P)+H(l,4)F°(q,P)= / k 3 dk / sin 2 Mf(l, 1), 

JQ Jo 

H(2,l)F°(q,P) + H(2,2)F°(q,P) + H(2,3)F^q,P)+H(2,4)F°(q,P) = / k 3 dk / sin 2 (2, 2) + K(2, 3)), 

/'OO /»7T 

1/(3, P)+P,(3,2)P 2 °(g,P) + /P(3, 3)^, P) + P (3, 4)P 4 °(g, P) = / fc 3 dfc / sin 2 (3, 2) + X(3, 3)), 

Jo Jo 

/■ OO /»7T 

i/(4,l)F{ } ( g ,P)+/P(4,2)F 2 (g,P) + i?(4,3)F3 1 ( (Z ,P) + P(4,4)F 4 ( g ,P)= / fc 3 dfc / sin 2 0ff(4,4), (12) 

Jo Jo 

the expressions of the H(i,j) and K(i,j) are cumbersome and neglected here. 

Here we will give some explanations for the expressions of H(i,j) . The H(i,j)'s are functions of the quark's 
Schwingcr-Dyson functions (SDF) A(q 2 + <f P 2 +£q-P), B(q 2 + £ 2 P 2 +£q-P), A(q 2 + (1 - £) 2 P 2 - (1 - £)q ■ P) and 
B{q 2 + (1 — £) 2 P 2 — (1 — 0? ' P)- The relative four-momentum g is a quantity in the Euclidean space-time while the 
center of mass four-momentum P is a quantity in the Minkowski space-time. The present theoretical techniques can 
not solve the SDE in the Minkowski space-time, we have to expand A and B in terms of Taylor series of g • P, for 
example, 

A(q 2 + eP 2 + &■■?) = A{q 2 + eP 2 ) + A{q 2 + £ 2 P 2 )'£g • P + • • • . 

The other problem is that we can not solve the SDE in the time-like region as the two point gluon Green's function 
can not be exactly inferred from the SU (3) color gauge theory even in the low energy space-like region. In practical 
manipulations, we can extrapolate the values of A and B from the space-like region smoothly to the time-like region 
with suitable polynomial functions. To avoid possible violation with confinement in sense of the appearance of pole 
masses g 2 = —m{q 2 ) in the time-like region, we must be care in choosing the polynomial functions [17]. 
Finally we write down the normalization condition for the BSW, 

/ (0<*^^ - (' - ^ + ^ + ™«- p ^at° P> i " 2P - (13) 



where x — liX + H- We can substitute the expressions of the BSWs and SDFs into the above equation and obtain the 
analytical result, however, the expressions arc cumbersome and neglected here. 

V. COUPLED RAINBOW SD EQUATION AND LADDER BS EQUATION AND THE DECAY 

CONSTANTS 

In this section, we explore the coupled equations of the rainbow SDE and ladder BSE for the pseudoscalar mesons. 
In solving those equations numerically, the simultaneous iterations converge quickly to an unique value independent 
of the choice of initial wavefunctions. The final results for the SDFs and BSWs are plotted as functions of the square 
momentum q 2 . 

The quark-gluon vertex can be dressed through the solutions of the Ward-Takahashi identity or Salvnov- Taylor 
identity and taken to be the Ball-Chiu vertex and Curtis-Pcnnington vertex [24,25]. Although it is possible to solve 
the SDE with the dressed vertex, our analytical results indicate that the expressions for the BSEs with the dressed 
vertex are cumbersome and not suitable for numerical iterations **. 

In order to demonstrate the confinement of quarks, we have to study the analyticity of SDFs for the u, d, s, 
c and b quarks, and prove that there no mass poles on the real timelike q 2 axial. It is necessary to perform an 
analytical continuation of the dressed quark propagators from the Euclidean space into the Minkowski space q^ — ► iqo. 
However, we have no knowledge about the singularity structure of quark propagators in the whole complex plane as 
our solutions are obtained in the Euclidean regions. We can take an alternative safe procedure, stay completely in 
the Euclidean space and take the Fourier transform with respect to the Euclidean time T for the scalar part (S s ) of 
the quark propagator [6,11,26], 



where the 3-vector part of q is set to zero. If S(q) has a mass pole at q 2 = —m 2 (q 2 ) in the real timelike region, the 
Fourier transformed S*(T) would fall off as e~ mT for large T or logS** = —mT. 

In our numerical calculations, for small T, the values of S* are positive and decrease rapidly to zero and beyond 
with the increase of T, which are compatible with the result from lattice simulations [27] ; for large T, the values 

"This observation is based on the authors's work in USTC. 




(14) 
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of 5* are negative, except occasionally a very small fraction positive values. We can express S* as |5*|e m7r , n is an 
odd integer, log 5** = log 15**1 + inn. If we neglect the imaginary part, we find that when the Euclidean time T is 
large, there indeed exists a crudely approximated (almost flat) linear function with about zero slope for all the u, d 
(in isospin symmetry limit, the u and d quarks are equal), s, c and b quarks with respect to the variable T, which 
is shown in Fig.l. Here the word 'crudely' should be understand in the linearly fitted sense, to be exact, there is 
no linear function. However, such crudely fitted linear functions are hard to acquire physical explanation and the 
negative values for S* indicate an explicit violation of the axiom of reflection positivity [28], in other words, the quarks 
are not physical observable i.e. confinement. 

From Fig. 2, we can see that for the u, d (in isospin symmetry limit) and s quarks, the SDFs are greatly renormal- 
ized at small momentum region and the curves are steep at about q 2 = lGeV 2 which indicates an explicit dynamical 
chiral symmetry breaking, while at large q 2 , they take asymptotic behavior. As for the c and b quarks, shown in 
Fig. 3, the current masses are very large, the renormalization is more tender, however, mass poles in the time-like 
region are also absent, which can be seen from Fig.l. At zero momentum, m u (0) = 688MeV, m d (0) — 688MeV, 
m a (0) = 882MeV, m c (0) = 1823MeV and m b (0) = 4960MeV r , which are compatible with the Euclidean constituent 
quark masses defined by m 2 (q 2 ) = q 2 in the literature. From the plotted BSWs (see Fig. 4 as an example), we can 
see that the BSWs for pseudoscalar mesons have the same type (Gaussian type) momentum dependence while the 
quantitative values are different from each other. The gaussian type BS wavefunctions which center around small 
momentum indicate that the bound states exist in the infrared region. Finally we obtain the values for the decay 
constants of those pseudoscalar mesons which are defined by 



,U = \27MeV\ f K = 157MeV; f D = 238MeV: f Ds = 2AlMeV: f B = 193MeV; f Bs = 195MeV, (16) 



which are compatible with the results from the experimental extractions, lattice simulations and QCD sum rules 



calculations, f„ = 130MeV(Exp), f K = 160MeV{Exp), f D = 300l^t^MeV r ( J Ba;p), f Ds = 285±19±40MeV{Exp), 
f D w 160 - 230MeV{Latt, sumrule), f Ds w 210 - 260MeV(Latt 1 sumrule), f B « 150 - 230MeV{Latt, sumrule) 



iUPy. = (0\qWaq\*(P)) 




d 4 k ' 



(15) 



here we use tt to represent the pseudoscalar mesons. 
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and Jb s ~ 190 — 260MeV(Latt, sumrule) [29-31]. The values from the lattice simulations and QCD sum rules 
calculations vary in a broad range, here we make an approximation. In calculation, the input parameters are N = 1.0A, 
V(0) = -11. OA, p = 5.0A, m u = m d = 6MeV, m s = ISOMeV, m c = 1250MeV, m b = 4700MeV, A = 200MeV, 
vj = l.6GeV and A = O.OAGeV 2 , the masses of the pseudoscalar mesons are taken as input parameters [21]. 

VI. CONCLUSION 

In this article, we investigate the under-structures of the pseudoscalar mesons (ir, K, D, D s , B and B s ) in the 
framework of the coupled rainbow SD equation and ladder BS equation with the confining effective potential (infrared 
modified flat bottom potential). Although the quark-gluon vertex can be dressed through the solutions of the Ward- 
Takahashi identity or Salvnov- Taylor identity and taken to be the Ball-Chiu vertex and Curtis-Pennington vertex, a 
consistently numerical manipulation is unpractical, we take bare approximation for the vertex T M = 7 M . After we solve 
the coupled rainbow SDE and ladder BSE numerically, we obtain the SDFs and BSWs for the pseudoscalar mesons. 
The SDFs for the u, d and s quarks are greatly renormalized at small momentum region and the curves are steep at 
about q 2 = IGeV 2 which indicates an explicitly dynamical chiral symmetry breaking. After we take the Euclidean 
time fourier transformation about the quark propagators, we can find that there are no mass poles in the time-like 
region, and obtain satisfactory results about confinement. As for the c and b quarks, the current masses are very 
large, the renormalization is more tender, however, mass poles in the time-like region are also absent. The BSWs for 
the pseudoscalar mesons have the same type (Gaussian type) momentum dependence while the quantitative values 
are different from each other. The gaussian type BS wavefunctions which center around small momentum indicate 
that the bound states exist in the infrared region. Our numerical results for the values of the decay constants of 
the pseudoscalar are compatible with the corresponding ones obtained from the experimental extractions and other 
theoretical calculations, such as lattice simulations and QCD sum rules. Once the satisfactory SDFs and BSWs for 
the pseudoscalar mesons are known, we can use them to investigate a lot of important quantities in the B meson 
decays, such &sB — -k,B — K,B — D,B — p former factors, Isgur-Wise functions, strong coupling constants, etc. 
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Fig. 3 SDFs of c and b quarks 
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Fig.4 BSVvs of n meson 



FIG. 1. Fig.l -Log\S*.(T)\. 
FIG. 2. Fig.2 SDFs of u and s quarks. 
FIG. 3. Fig.3 SDFs of c and b quarks. 
FIG. 4. Fig.4 BSWs of tt meson. 
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